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Using a recently proposed classification scheme for phase transitions in finite systems 
[Phys. Rev. Lett. 84, 3511 (2000)] we show that within the statistical standard model of nuclear mul- 
tifragmentation the predicted phase transition is of first order. 
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I. INTRODUCTION 



of a first order phase transition are presented in Sec. IV 
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The phenomenon of heavy-ion collision induced nuclear 
multifragmentation (NMF) has been extensively studied 
over the past decade. Experimentally it is found that 
nuclei with excitation energies of E* /A = 5-10 MeV, 
i.e. in the region of the binding energy, expand forming 
a number of intermediate mass fragments M . Using sub- 
tle and advanced experimental techniques from a num- 
ber of collision experiments with different nuclei and at 
different collision energies an experimental caloric curve 
has been derived p. p|, which supports the interpreta- 
tion that NMF is the small particle number onset of the 
nuclear liquid-gas phase transition. Theoretically NMF 
has been described by a number of different approaches 
ranging from simple percolation models [0, pi, dynamical 
models with different levels of sophistication, e.g. quan- 
tum molecular dynamics of different flavors [|[ 0, ||, ^, |I^ 
and fermion molecular dynamics pd|, a nd a large vari- 
ety of statistical models |l|, ^, |14|, [15| |l|, |l^, ^, ^. 
Among the statistical models two mainstream-theories, 
namely the statistical multifragmentation model (SMM) 
| jl2[ using the heat bath or canonical ensemble and the 
micro-canonical multifragmentation model (MMM) |13|| 
using the constant energy or micro-canonical ensemble 
should be mentioned. For different ensembles the shape 
of the caloric curve may differ significantly. However, 
because all ensembles are uniquely connected by simple 
integral transforms the qualitative features of the phase 
transition should be the same in all ensembles. We there- 
fore omit at this point the (surely) interesting question 
which ensemble is more appropriate for the description 
of nuclear multifragmentation. 

In this paper we use the canonical statistical multi- 
fragmentation model based on the works of Mekjian et 
0, |1|, 



|2l| , p3, E3| to analyze the order of 
the multifragmentation phase transition using a recently 
proposed classification scheme for phase transitions in fi- 
The classification scheme is based 
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nite systems 

on an analysis of the distribution of zeros of the canoni- 
cal partition function in the complex plane and is in the 
thermodynamic limit in complete accordance with the 
Ehrenfest classification scheme. This scheme is shortly 
described in Sec. II. In Sec. [II we briefly review the 
statistical model used in this paper. The results of our 
analysis which clearly show that NMF has the signature 



II. CLASSIFICATION SCHEME 

The classification scheme used below and described in 
far more detail in |gj, [2^ is based on the analytic continu- 
ation of the inverse temperature /3 into the complex plane 
f3 —> B = P + ir. For macroscopic systems an equivalent 
scheme was developed by Grossmann et al. [E6| . All ther- 
modynamic information about a system can be extracted 
from the distribution of complex zeros of the canonical 
partition function Z{B) ^ J dE n{E) exp{-BE). 

To avoid some difficulties with the high tempera- 
ture limit we write the partition function as a product 
Z{B) — Z\i^{B)Zint{B), exphcitly separating the high 
temperature limit (T —>■ cxd) of the partition function 
Ziim{B) imposing limT^oo Zint{B) = 1. The partition 
function is an integral function with complex conjugate 
zeros Bk = B*_f, = Pk + iTk {k 6 N) and can be rewritten 
as 



Z{B) = Zii,„(S) Zint(O) exp(SaBlnZint(0)) 
B\ /^ B\ f B 
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All thermodynamic properties like the internal energy 
or specific heat are derivatives of the logarithm of the 
canonical partition function. Thus, the zeros of Z{B) are 
poles of these quantities. In the thermodynamic limit 
different regions of holomorphy separate different phases 
by dense lines of zeros. In finite systems the zeros do 
not become dense on lines which leads to a less sharp 
separation of different phases. We interpret the zeros as 
boundary posts between two phases. 

The distribution of zeros close to the real axis can ap- 
proximately be described by three parameters, where two 
of them refiect the order of the phase transition and the 
third merely the size of the system. We assume that the 
zeros lie on straight lines with a discrete density of zeros 
given by 
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with k = 2,3,4, •••, and approximate for small r the 
density of zeros by a simple power law </'(t) ~ t". Con- 
sidering only the first three zeros the exponent a can be 
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FIG. 1: Schematic illustration of the zeros in the complex 
temperature plane. 



estimated as 
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The second parameter describing the distribution of ze- 
ros is given by 7 = tan v ~ ((32 — Pi) / {t2 — ti) where v is 
the crossing angle of the line of zeros with the real axis 
(see Fig. ||) . The discreteness of the system is reflected 
in the imaginary part ri of the zero closest to the real 
axis. For macroscopic systems we always have ti —> 0. 
In this case the parameters a and 7 match those defined 
by Grossmann et al. [^ , who showed the connection be- 
tween certain values of a and 7 and different types of 
phase transitions. They claimed that a first order tran- 
sition occurs for a = 7 = 0, for < a < 1 with 7 = 
or 7 7^ a second order transition, and for a > 1 a 
higher order transition. In the thermodynamic limit a is 
restricted to values greater than or equal to zero because 
otherwise this would cause a diverging internal energy. 
For finite systems values of a less than zero are possible. 
In our classification scheme phase transitions in finite 
systems are of first order if a < 0, while the definitions 
for higher order transitions coincide with those given by 
Grossmann. Additionally we define ri as an unambigu- 
ous parameter for the discreteness of the system. Since 
in small systems no thermodynamic properties diverge, 
the specification of the critical temperature is difficult 
and several definitions are possible which coincide only 
in the thermodynamic limit. We define the real part of 
the first zero /3i as the critical temperature. Another 
possible definition would be the crossing point of the line 
of zeros with the real temperature axis. 



of the nuclear force. The properties of the systems at 
this so-called freeze-out density, which is assumed to be 
0.1-0.3 of the normal nuclear density po = 0.15 fm~'^, are 
than carried by all fragments until experimental detec- 
tion. 

The probability for a fragmentation with fragmenta- 
tion vector n — (m, 712, . . .) {nk is the number of clusters 
of size k and ^ kn^ = A), reads 



w{n) - n 
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where x 



K, 



M 



\3/2 



is the partition function for 



the ideal monoatomic gas of mass M in the accessible 
volume T4c. y = exp (^ay{B)B + jfxk) ^^ related to the 
cluster binding and internal excitations. The cluster size 
dependence of the probability is given by a^ = fc^/^. 
We assume an accessible volume three times larger than 
the normal nucleus volume Vq. The canonical partition 
function is then given as 



\^ ^("^)^m^,A~m 



m=0 



(5) 



where m — ^^, Uk is the multiplicity, i.e. the total num- 
ber of clusters, and Z^' is independent of the thermo- 
dynamic variables and given by the simple recursion 



1 ^ 1 



(6) 



with Z)^ = 1/aA- 

Following Ref. |Q the binding energy for a composite 
of k nucleons is given by 



a„(S) = Wofc-cr(S)fc2/3, 



(7) 



with a{B) = ffo {{Tc - IIB'^)/{T^ + l/B'^)f'^ being the 
surface tension (Wo = 16MeV, ctq = 18MeV, and T^ = 
18MeV). 

Within this framework different moments of the mul- 
tiplicity are calculated by 



1 ^ 
i^m') = ^ m'Z 



Za 



M^m^A-™_ (8) 



III. STATISTICAL MULTIFRAGMENTATION 
MODEL 

In the statistical models for NMF [|l|, |l^, ||, |2| it is 
assumed that the fragmentation process can be described 
in the following manner. After initial excitation the sys- 
tem expands and fragments into a number of clusters. 
During the expansion a thermalization of all clusters is 
assumed. The thermalization abruptly stops when the 
mean separation between the clusters exceeds the range 



The variance of the multiplicity is given by Atti = {w?) — 
(m)'^. It is straightforward to express thermodynamic 
quantities like the internal energy or the specific heat in 
terms of different moments of the multiplicity |17 . 



IV. RESULTS 

Due to the unfavorable scaling behavior of the parti- 
tion function it is very arduous to find its zeros directly 
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FIG. 2: Contour plots of the variance of the multiphcity in the complex temperature plane for (a) A=25, (b) A=75, (c) A=150, 
and (d) A=225 nucleons. The dark spots indicate the locations of the zeros of Z{B). 



with numerical methods. The zeros of Z are the poles 
of almost all thermodynamic expectation values, e.g. the 
internal energy, the specific heat and so on. In this pa- 
per we utilize the absolute value of the variance of the 
multiplicity | (Am) | to locate the zeros of Z in the com- 
plex temperature plane. Fig. shows a contour plot of 
{Am)\/A'^ in the complex temperature plane for nucleon- 
numbers A = 25, 75, 150, and 225 (note the different scal- 
ing of the r-axis.). This figure directly shows that multi- 
fragmentation is a first order phase transition. The zeros 
lie on straight lines perpendicular to the real tempera- 
ture axis with constant spacings between them. Even at 
particle numbers as small as 25 the nature of the phase 
transition can be clearly identified from the lines of zeros 
of the partition function. In contrast, it is very diffi- 
cult to extract the order of the phase transition from the 
canonical caloric curve (internal energy Uint per particle 
number A vs. temperature T) shown in Fig. 0. 

The recursion formula (0) has a very similar structure 
to that for Bose-Einstein systems [g5[ ^, |8|. This is 
because both systems have an intimate relationship to 
the permutation group |23]. It is therefore remarkable, 
that the multifragmentation phase transition for nuclear 
matter is of first order while Bose-Einstein condensation 
is a second or third order phase transition depending on 
the single particle density of states. 

A detailed analysis of the distribution of zeros reveals 
that for all nucleon numbers {A > 25) the most impor- 
tant classification parameter a is almost equal to zero 
(see Fig. ||(a)). The slight deviations for A=50 and 
A=150 are most likely due to uncertainties within our 
numerical zero detection method. The second parameter 
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FIG. 3: Canonical caloric curve for A=25, 100, and A=225. 
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FIG. 4: Nucleon number dependent value of the classification 
parameters (a) a, (b) r for an accessible volume 14c = 3Vo and 
(c) of the critical temperature Pi for three different accessible 
volumes. 



7 can directly be extracted from Fig. 0. Even for small 
nucleon numbers the distribution of zeros is perpendicu- 
lar to the real axis, corresponding to 7 = 0. The finite 
size of the system is connected to the third classification 
parameter ri. Fig. H(b) shows the yl-dependence of ti. 
With increasing nucleon number ti decreases. We ap- 
proximated the scaling behavior as ti ~ ^-0-96^ [ g ^^ jg 
simply volume dependent and the phase transition would 
approach a true first order phase transition in the Ehren- 
fest sense with 1/A. The critical temperature T = l//3i 
(Fig. ^c)) decreases and adopts at A « 200 the exper- 
imentally found value of about 5 MeV. Of course, for 
larger system sizes the simple equation (0) for the bind- 
ing energy used in this study should be refined and aug- 
mented by additional terms. Additional calculations for 
different particle densities reveal that there are no quali- 



tative changes to the first order nature of the phase tran- 
sition. Fig. ^(c) also shows calculations for accessible 
volumes which are twice, three, and five times as large as 
the normal nucleus volume Vq. There are only quantita- 
tive changes in the critical temperature. 

It should be noted that one can expect that refine- 
ments of the model, e.g. the introduction of an additional 
Coulomb term and other values of the particle density, 
may influence the details of the distribution of zeros. Es- 
pecially tuning the particle density to the expected su- 
percritical phase is an interesting challenge, worth to be 
pursued in the future. 



order phase transition, which is in complete agreement 
with recent works based on micro-canonical and grand- 
canonical calculations [^, Q. Of course, the first order 
type of the NMF phase transition is not yet sufficient 
to intcrprete NMF as the onset of the nuclear liquid-gas 
transition. However, this seems to be a necessary condi- 
tion. The distribution of zeros of the partition function 
in the complex temperature plane proofs as a reliable tool 
for the classification of phase transitions, adding a clear, 
detailed, and unambiguous view on the thermodynamic 
properties of small systems. 



CONCLUSION 
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